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Abstract 

The covariant kinetic approach for the radiative plasma, a mixture 
of a relativistic moving gas plus radiation quanta (photons, neutrinos, 
or gravitons) is generalized to D spatial dimensions. The operational 
and physical meaning of Eckart's temperature is reexamined and the 
D-dimensional expressions for the transport coefficients (heat conduc- 
tion, bulk and shear viscosity) are explicitly evaluated to first order 
in the mean free time of the radiation quanta. Weinberg's conclusion 
that the mixture behaves like a relativistic imperfect simple fluid (in 
Eckart's formulation) depends neither on the number of spatial dimen- 
sions nor on the details of the collisional term. The case of Thomson 
scaterring is studied in detail, and some consequences for higher di- 
mensional cosmologies are also discussed. 
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I. INTRODUCTION 



In many astrophysical and cosmological problems, the matter-radiation content 
is modeled by a so-called radiative plasma, a two component fluid consisting of 
some material medium (perfect fluid, relativistic gas, ionized hydrogen atoms, etc.) 
plus radiat ion quanta (photons, neutrinos, or gravitons). In the simplest situation, 
the material medium is assumed to be locally in thermal equilibrium with itself 
(very short mean free times and mean free paths), whereas the massless component 
has a finite mean free time r . However, even in this case, an idealized perfect 
fluid description of this mixture is inappropriate, and a deeper insight is provided 
by a kinetic theory approach. In particular, since the radiation quanta are out of 
equilibrium with the material component, all transport properties of this system, to 
first order of approximation, depend exclusively on the mean free time of radiation. 

The first self-consistent application of the relativistic kinetic theory to this sys- 
tem is due to Thomas 0, who computed the correct values of the heat and shear 
viscosity transport coefficients. Later on, by considering the concept of local equi- 
librium temperature, Weinberg obtained the bulk viscosity coefficient, thereby 
showing that such a mixture behaves like an imperfect relativistic simple fluid in the 
framework of the hydrodynamical formulation for dissipative processes developed 
by Eckart 0. Masaki ||, recovered the results of Thomas using a modern rela- 
tivistic approach (manifestly covariant). This work was considerably simplified by 
Straumann [||] , who also shown that if the Thomson scattering dominates, the bulk 
viscosity coefficient goes to zero. More recently, by considering the 9-moment Grad 
approximation method to solve the Boltzman equation, Schweizer calculated the 
transient transport coefficients appearing in the linear, causal, hydrodynamical for- 
mulation proposed by Israel and Stewart . The physical and operational meaning 
of local equilibrium temperature (Eckart 's temperature) has also been discussed by 
several authors 0, [|], 0. 
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On the other hand, since multidimensional models are presently believed to be 
the natural framework to unify all fundamental interactions, it is very common nowa- 
days to study all phenomena and the underlying physics in an arbitrary number of 
dimensions. Historically, arbitrary dimensionality (including fractal dimension) has 
been widely considered in the context of Newtonian gravitational theory and gen- 
eral relativity, quantum field theory, cosmology, phase transitions, renormalization 
group, quantum mechanics, lattice models and random walks pl j. In domain of 



cosmology, for instance, several compactification mechanisms for extra dimensions, 
presumably taking place in the very early universe, have been suggested to explain 



the now observed 4-dimensional structure of spacetime [|T^ . 

Given the lack of knowledge about the initial state of the universe, the sim- 
plest dynamic approach to dimensional reduction is to construct a multidimensional 
anisotropic cosmological model based on some fundamental theory (superstrings, 
Kaluza-Klein or Einstein's gravitational theory) where the extra spatial dimensions 
are collapsing while the "physical dimensions" are expanding, and eventually, leading 
to the present day universe. Isotropization and cosmological dimensional reduction 
has also been considered in the framework of nonequilibrium higher dimensional 
cosmologies. Concrete models have been proposed in the case of bulk and shear 
viscosities leading to the conclusion that viscosity may work efficiently to reduce the 



number of spatial dimensions. [|Tj]- fl6| . As a matter of fact, this is a new version 
of the earlier program of chaotic cosmology, for which the ordinary dissipative pro- 
cesses (viscosity, heat flow, ...etc.) were believed to be the mechanisms damping out 
any anisotropy and inhomogeneity existing in the primeval plasma [[17]]. However, 
even considering that the kinetic approach is more fundamental than its related 
thermodynamics, the analysis mentioned above are usually based only on macro- 
scopic equations. In particular, to the best of our knowledge, the explicit form of the 
transport coefficients considered in such works have not been derived from kinetic 
theory in D spatial dimensions. This one of the main aims of this article. 



It is also worth noticing that the definition of a radiative plasma, in princi- 
ple, is not restricted to the class of systems mentioned above. Under certain con- 
ditions, more exotic system as, for instance, a mixture of massive and massless 
modes in superstrings models could, formally, be interpreted as a radiative plasma 
in D dimensions. The reason for this possibility may be easily understood. In the 
collisional limit, two different temperatures are naturally defined in an radiative 
plasma, namely: the matter temperature(T m ) and the effective radiation kinetic 
temperature (T r ). However, to first order in the mean free time, the latter is differ- 
ent from the former only if the exp ansion rate is different from zero. The expansion 
works to continuously pulling radiation and matter out of thermal equilibrium, be- 
cause the components have different cooling rates. Basically, this is the mechanism 
accounting for the presence of the bulk viscosity in such a mixture as well as in a 
relativistic simple gas ||. On the other hand, even in a thermodynamic setting, a 
gas of hot strings(superstrings, or Hagedorn's fireballs) may also be described as a 
two temperature system when we have a mixture of heavy strings with a gas of light 
strings(see, for instance, [0). 

From the above considerations, it seems interesting to analyze the transport 
properties of a D-dimensional radiative plasma. In this work, by extending the 
covariant 4-dimensional kinetic approach developed by Straumann ||, we compute 
the transport coefficients for such a system. Hopefully, the expressions derived here 
may be useful to study the effects of dissipative processes in higher dimensional 
cosmologies as well as in modeling a mixture of massive and massless modes in the 
framework of superstring theory. 

The paper is structured as follows: In the next section, the energy-momentum 
tensor for a D-dimensional radiative plasma is computed exactly by solving the 
D-dimensional Boltzmann transport equation to first order in the mean free time 
of the radiation quanta. In section 3 we discuss the physical meaning of Eckart's 
temperature and obtain the transport coefficients of heat conduction, bulk and shear 
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viscosities, and in section 4, the photon and total entropy production are discussed. 
In section 5, a specific model, namely, Thomson scattering, is studied. We conclude, 
in section 6, with a discussion of our main results for higher dimensional cosmology, 
and in the following two appendix we proof the basic mathematical properties used 
in the article. In our units the physical constants h = kg = c = I, where ks is 
Boltzmann's constant and the D+l-dimensional signature of the spacetime metric 
is (+,-,-,-,...,-). 



II. D-DIMENSIONAL KINETIC APPROACH 

The basic object for computing the energy-momentum tensor(EMT) and 
other macroscopic quantities describing a D-dimensional radiative plasma is the 
radiation distribution function F(k,x). As usual, it is defined in such a way that 
F(k, x)d D xd D k gives the number of quanta in the D-dimensional volume element 
located at x, and whose D-momentum k lie within d D k. Now, by assuming that 
the radiation quanta (photons for the sake of simplicity P2| ) are out but close to 
equilibrium, we may expand F(k,x), to first order in mean free time, as 

F = F (o) + F (i) ) (1) 

where F^ is the equilibrium distribution in the D-dimensional space 

F<°> = — - (2) 

(27r)" e ¥-l ' [) 

and the interaction part F^ satisfies \F^\ « F^ '. In the above expressions, is 
the normalized four-velocity field {U^U^ = 1) of the material medium and T = T(x) 
its local temperature which, in equilibrium, coincides with the uniform temperature 
T M of the material medium. 

The distribution function satisfies the Boltzmann transport equation which to 
first order may be written as 
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dF {0) 



L[F(% 



(3) 



dx^ 



where L is a linear operator which depends on the specific model (interaction of 
radiation and material fluid). Following Strauman ||, the solution of the above 
equation will be discussed, as far as possible, in rather general terms. An application 
of these results for a specific model will be presented in section 5. 

As is widely known, the above Eq.(|3]) cannot, in general, be exactly solved so that 
some approximation scheme need to be implemented. One which is often considered 
is provided by the Grad's relativistic 9-moments method applied in Refs. [[||| to the 
4-dimensional case. In the present D-dimensional framework it may be generalized 
as is explained below. 

Let G be the D-dimensional Lorentz group of the Minkowski-like space-time M 
and G x the little group associated to a timelike vector U(x). This subgroup of 
G leaves U(x) invariant. As a consequence G x is isomorphic to SO(D), thereby 
implying that F(k, x) is a function of u = k^U 11 , n M and x where 



Note that n^U^ = and n ll n fl = — 1. As we shall see next, the key ingredient to 
make the integrals on the hypersphere S D _ 1 (space-group where G x acts upon) easier 
to compute comes from the fact that u = k^U^ is obviously invariant under G x at 
every point x. 

The calculational basis of Grad's method applied to radiative problems is to 
expand the perturbative term F^(w, x) with respect to n M into irreducible poly- 
nomials under the action of the little group G x 



F^\uj,n^x)=A(uj,x) + B^uo,x)n^ + C^(uo,x)(n^n u + -h^) + ... , (5) 




(4) 



where 



V 




(6) 
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is the projector onto the hyperplane normal to £/ M , and jf v is the Minkowski matrix. 
The former object also satisfies the identities: h%U" = and h^n" = n M . Naturally, 
the vectorial and tensorial parts in the decomposition (5) are defined in the rest space 
normal to £/ M i.e., h%B v = B^ and h°h*C\ a = C^ u . In addition, since n^n u + jjh£ is 
symmetric and traceless, without loss of generality, one may assume that the same 
holds for C^ v , that is: 

Cj//i — , = . (7) 

On the other hand, since L behaves like a scalar under the action of the little 
group G x , this means that it operates in the irreducible subspaces spanned by the 
irreducible polynomials in (5) as a multiple of the unit operator. It thus follows that 

L[F&] = -u[k A + k ± BX + KtC^inW + ^/O + ...] , (8) 

where = r" 1 (i = 0, 1, 2), is the inverse of the mean free time for the process in 
consideration (bulk viscosity, heat flux, shear viscosity). As usual, we assume that 
these quantities are functions solely of cu and x. 

Using the transport equation (3) we can express the coefficients A, B^ and 
in terms of £/ M , T, and To do that we introduce a measure dVLu on the (D — 1)- 
dimensional hyper surf ace 

S D1 = {k/ k 2 = 0, k° > 0, k.U = const} , (9) 

for which dVLjj is the unique G x -invariant measure normalized as 

r 1-n Dj2 

L M u = - m = s D _ 1 , (10) 

Jb D-l 1 I 2 I 

where T is the gamma function. 

This normalization is easily understood in the comoving frame. Its nothing 
more than the "area" of unit (D — l)-dimensional sphere. In the above formula we 
have denoted the area of S D1 with this same symbol, but this willnot cause any 
misunderstanding to the reader. 
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Now, with respect to this measure, the irreducible polynomials are normalized 
by (see Appendix 1) 



S D-1 D — l D 

J— I {n»n u + ±-W v ){rfn p + hf p )d£lv = 

^D-l JS D-1 U U 

— — AhTh"" + h^h vp + h^h va ) - -^h^h ap , 

D f D I 2 ) 



;n) 



(12) 



Also we have the very useful result 



/ n»n v . . . rfdVtu = , 
Js„ , 



(13) 



for an odd number of components of the vector n. From a straightforward 
(although unduly long) calculation we can get the coefficients A, B^, C^ u in (5) 
by evaluating the moments of the Boltzmann equation (3), and using the relations 
(11)-(13). The results are: 



1 UJ , I UJ 



T 

- 1 - 71. 



—U p d T + — 



(14) 



B„ 



1 UJ 
K l T M 



UJ 



T, 



-h^d v T M - U»d u U, 



(15) 



C -_J_^. 



where 4>{w~) = F°(uj, x), and a prime denotes its derivative. 

M 

By definition, the total energy-momentum tensor T pv of the mixture (matter 
plus radiation) is given by 

, d D k 



rpfjLP rpflU _|_ rpfMU 



(17) 



where is the EMT of the material component, which from now on, will be 
described by a perfect fluid 



T& = p M U"ir - , (18) 

and p M and p M are, respectively, the energy density and pressure. Tj^ is the radia- 



tion EMT for the equilibrium state, which may be readily evaluated to be 

Tg = J k^F ^- = a D T^ (W - , (19) 

or equivalently, in the standard form 

T]? = p^Wir - p H hr , (20) 



where 



a D T M + \ Pr = JjPR 



and 



4 7T D / 2 



-n- . Dr(D S(D + i)((D + i) , (21) 

is the D-dimensional radiation constant, and ( denotes the Riemann zeta function. 

The last term (the integral) in Eq.(17) is the only out of equilibrium contribution 
(due to interaction between matter and radiation), and can properly be interpreted 
as a dissipative term to the equilibrium energy-momentum tensor. The coefficients 
A, Bp, C^v give rise to contributions of different tensorial ranks, which in the ther- 
modynamic theory are usually named, respectively, as bulk viscosity, heat flux, and 
shear tensor. Substituting (14)-(16) in the last integral of (17), we get the terms: 

a) Scalar (bulk viscosity) F^/ = A 



k° So D \T M A M D 



b) Vectorial (heat flow) F$ = B^n^ 



■ x 



(22) 



J k»k»Fff^Jl = (^^) ^k(h^U" + h vX U»){d x T M - T M U°U x , a ) (23) 



c) Tensorial (shear viscosity) F^ h ' = C fiu (n^n u + -kh^) 

I kVF ^ - (J WrSr mu '-' ~ v >* ~ 2 d k ' v * ] • (24) 

In the above equations, Ri (i = 0, 1, 2) are the D-dimensional Rosseland means for 
Ki. In general, for an arbitrary quantity we define its D-dimensional Rosseland 
mean by 

Ki ft* U>l>+W^)du> ' 1 ) 

M 

III. ECKART TEMPERATURE AND THE TRANSPORT COEFFICIENTS 

In the phenomenological theory of Eckart, the EMT of a relativistic dissipative 
simple fluid reads 

= p(T, n)WU v - p(T, n)h^ u + AT"" , (26) 

where the last term is the macroscopic description of the irreversible processes. Its 
D-dimensional canonical form is given by (the factor 1/D appears explicitly because 
the shear viscosity stress is traceless) 

AT"" = ( D W U U X X + x D {W X U u + h vX U»)(T x - TU a U x ,a) 
+VD h^h^(U a , p + U p , (T -^g ap U x x ) , (27) 

where ( D , Xm an d Vd are > respectively, the bulk viscosity, heat conducting, and 
the shear viscosity transport coefficients. In principle, the above expression should 
be justified by kinetic contributions like (22)-(24). It is worth noticing that all 
thermodynamic quantities (explicitly or implicitly) present in (26) are measured in 
the rest frame of U a . However, for a radiative plasma there are, in the collisional 
limit, at least three concepts of temperature ||. Since the physical features relating 
these temperatures concepts are independent (to first order in Rq x ) of the shear and 
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heat conducting effects, for a moment, we will restrict our discussion to the isotropic 
case for which only bulk viscosity is relevant. To do that, let us first introduce (from 
the bulk viscosity contribution (22)) the following auxiliary quantity: 

II M = —{D + l)a D T^- (^U x d x T M + ±U X X ) . (28) 

In terms of I1 M , the isotropic part of the total EMT (17) can be written as (see (22)) 

T"* = [ PM (T M ,n) + a D T° +1 + IIjEW - \p M (T M ,n) 

+^ T^ +1 + ^nj/^ . (29) 



K,T [^rU x d x T M + ^U\. 



(30) 



Note that by defining the effective radiation kinetic temperature: 

Tr = T M 

the above EMT assume the form 

TV = (p M (T M , n) + a D Tj^ +1 ) - (p M {T M , n) + ^a D T^ hT , (31) 

which has only two standard terms, which are characteristic of the perfect fluid 
description. In particular, the radiation quanta behaves as if they were in equilibrium 
at temperature Tr (compare (19)). However, the difference in temperature between 
the material and radiative components give rise to an irreversible heat transference 
between them. This is the physical mechanism accounting for the bulk viscosity 
process in this mixture. In order to extract its analytical expression, we recall that 
the Eckart temperature T is defined as being the local equilibrium temperature, 
which is implicitly fixed by (see Ref. 0) 

U,U V T^ = p(T, n) = p M (T M ,n) + a D T° +1 + U M = p M (T, n) + a D T D+1 . (32) 

Since the equilibrium temperature T e [T M ,T R ], whose length \T R — T M \ is of first 
order in Rq 1 , we can expand p(T,n) as a power series 

p(T,n) = p(T M ,n) + (J^j (T-TJ + --- , (33) 
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where 

p( T m , n ) = Pm(T m ,ti) + a D T^ +1 

Inserting (34) into (33) and comparing with (32) we get 

/ dp x 



(34) 



t = t + — n 



(35) 



Analogously, the total pressure p(T, n) can be expanded, to first order, as 
p (T,n)=p(T M ,n) + (^j (T-TJ + --- , 

or still, inserting p(T M ,n) from (31), and using (35) 

p{T,n) = p M (T M ,n) + ±a D T A D I +1 + n M . 



(36) 



(37) 



With the help of (28), (32), and (37), the total energy-momentum tensor (29) takes 
the form 

1 2 



T» v = p(T, n)U"U u - ^ p(T, n) 



1 I dp 



U\ } h» v (38) 



where we have used the thermodynamic relation 



T 



(39) 



and the fact that, to first order, T M can be replaced by T when multiplied by k . 
The above EMT has exactly the Eckart's canonical isotropic form, namely: 



= p(T, n)U»U u - [p(T, n) - (J)^ , 
where 9 = U\ is the expansion parameter, and from (38) and (40) we read 

2 



(40) 



C D = — a D T 



D+1 



1 

D 



dp 
Pp, 



(41) 



To obtain the remaining coefficients it is necessary only to observe that the vec- 
torial and tensorial contributions (23) and (24) are not modified (to the first order) 
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when rewritten in terms of the Eckart temperature. Hence, a direct comparison 
with the canonical form (27) give us 

X = D + la ° TD (42) 



and 

(D + l)a D T D+l 



(43) 



D(D + 2)R 2 ' 

Expressions (41), (42), and (43) are the transport coefficients for a D-dimensional 
radiative plasma. In particular, for D = 3, the standard results are recovered [0,0. 
From (41) we see that structureless point particles in D-dimensions have negligible 
bulk viscosity in the extreme relativistic regime (p = -kp), as should be expected 
in physical grounds. This is also a consequence of the fact that the D-dimensional 
equilibrium distribution of photons is preserved under expansion, or equivalently, 
that the ratio — is an adiabatic invariant in the sense of Ehrenfest, regardless of the 
number of spatial dimensions. Hence, when the universe is radiation dominated, 
there is no radiative bulk viscosity in the D-dimensional Friedmann-Robertson- 
Walker(FRW) spacetime. This generalize the widely known 3-dimensional result 



IV. ENTROPY PRODUCTION FOR A D-DIMENSIONAL RADIATIVE 

PLASMA 

In Eckart 's frame, the current of entropy and the local entropy production for a 
D-dimensional relativistic simple fluid are given by 

S» = s(T,n)U» + ^ , (44) 

and 

n 2 a n W v 

W* = ^ + ^ + ~T^ ' (45) 

Cd T X d T 2 2 Vd T 
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where s is the entropy density, and II, g M , and IP" are the classical dissipative 
stresses (bulk viscosity, heat flow and shear viscosity): 

n = c B , (46) 

q» = X D h»\T x -TU°U^) , (47) 

IP" = V D h^h^(U a , p + U Pia -^h ap U x x ) . (48) 

In principle, the Eckart description for a D-dimensional radiative plasma will be 
considered complete only if the kinetic equations lead to (44) and (45) in Eckart 's 
frame. In particular, one may argue that the above defined entropy density should 
be introduced, in a consistent way, by the same kind of algorithm used to fix the 
Eckart temperature. Indeed, this is exactly what happens. To show that let us 
consider the kinetic expression for the entro py flux || 

S£ = -j^f^[FlnF-(l + F)Hl + F)} , (49) 

where from (2), (14), (15), and (16) the distribution function F given in (5) may be 
written as 

f = * (f ) + Ko -f 0' (f ) n + « r l f a ( f ) 

\ M J M \ M / M \ M J 

KJ ' " '^^n^v + V) . (so) 



2 T M \T M ) ^ D 



where 



7 L / "^//', • j-r*; m 



% = ^rKd v T M - U»d u U, (52) 



IV = ^c/,, A + Ku^x - -K„u x x (53) 
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Inserting the above expression into (49) a straightforward, although lengthy, calcu- 
lation yields 

where has been indexed to recall that it is in the matter temperature. However, 
since q^ M is already of first order in R^ 1 , this means that -if- = C, and using the 

M 

definition of Tr given by (30), the above equation becomes 

S£ = ^a D T°U>> + £ . (55) 

Note that the radiation entropy density is U^S^ = — ^ should be expected 

in physical grounds. In particular, for D = 3, one has s# ~ T^, which, in the absence 
of heat flow, corresponds to the usual radiation fluid with p R ~ T^. These results 
are in agreement with our discussion in the earlier section. The quanta behave (in 
the isotropic case) as a radiation perfect fluid at temperature T R regardless of the 
number of dimensions. 

Now, as a self-consistency check of the algorithm used to introduce the temper- 
ature of Eckart, we need only to show that the total entropy density of the mixture 
reduces to the equilibrium expression at the temperature T. 

The total entropy of matter plus radiation is given by = S 1 ^ + S R where S R 
is given by (55) and = s M (T M ,n)£/ M is the entropy flux of matter. Since the 
material component is in equilibrium at temperature Tm, the total entropy density 
can be written as (see (54)) 

UpS* = s M (T„,n) + !^a D T° + ^ , (56) 

M 

where s M (T M , n) is the entropy density of the material medium. By definition of 
Eckart's frame, the right hand side of the above equation should be s(T, n). In fact, 
expanding s(T,n) in power series 

s(T,n) = s(T M ,n)+(^j (T — T M ) + o(Rq 2 ) , (57) 
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and using (35), the above expression (to first order) becomes 



s(T M ,n) = s(T,n)-(^) IT M , (58) 



where 



s (T M ,n) = s M (T M ,n) + I ^a D T° . (59) 

Inserting (58) into (56), we can rewrite the entropy density as 

= s(T,n) - ^ n M + ^ , (60) 

where we have used that Sr 4 - = ^M- to first order in f = Rq 1 . Finally, by considering 
the standard relation (||) = ^, we obtain 

U^ = s(T,n) . (61) 

Thus, in Eckart's temperature, the total entropy four-vector of a radiative plasma, 
assume the same form usually applied for a relativistic simple fluid (see (44)). 

Now, we proceed to calculate the photon entropy production, namely, d^S^. 
From (49) we have 

d^ R = -j^J d ^k^[F\nF-{l + F)Hl + F)} . (62) 

Using the Boltzmann equation (3) a straightforward calculation leads to 



where = 1 + F^. Using (50) these integrals are readily evaluated 
d^ R — (D + l)a D T D 



~ n 2 qjf i ruir^ 
n + — + ^f- + 



(64) 



R DRi 2D(D + 2) R 2 

where the linear term came from the first integral in (63). As it will be seen, if we 
include the contribution from the material medium to the total entropy production, 
we render the total linear term zero as a consequence of the conservation of the total 
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energy-momentum tensor (such approach has been applied by Schweizer || for the 
3-dimensional case). In order to do this we first introduce the notation 

% = r^ = r/3 v ■ (65) 

and using the Boltzmann equation (3), (65), and (19), the first integral in (63) can 
be evaluated leading to 



d^ R = (3 l/ T^ + QP , 



(66) 



where QP stands for the second integral in (63) which gives rise to the quadratic 
part in (64). Similarly, for the material medium one has 



f) CM _ ft rpflU 



(67) 



Adding (66) and (67), and using the conservation of the total energy-momentum 



tensor, T^ v = T R + T^, we get 



= QP=(D + l)a D T 



D 



R Dk x 2D(D + 2) R 2 
With a bit more algebra one can prove that 

n 2 1 n 2 

(D + l)a n T D — = ±thL , 



(68) 



(69) 



D 



«i X D T 2 



(70) 



(p + i)a D T D ~ ~ ryip- 

2D{D + 2)R 2 " u " 2 Vd T 



(71) 



where U^ v = rj D Ii^ v . Finally, using (69)-(71) we can write (68) in the expected final 
form 



in 2 i i 
duS" = ^-- + —^q^ + — ^n^n^ 



(72) 
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V. MODEL WITH THOMSON SCATTERING 



In D dimensions, the transport equation for the case of a nondegenerated material 
medium and temperature not too high, may be written in the form (see, for example, 
references [5], [6], and [14]) 



fcA»0 F = -uNa a (l - e^ /T )(F - F (0) ) 



—uNa, 



F — 



p(n, n')F(u, n')dVt v 



(73) 



In this expression, the quantities (y a {uj) and a s {uj) appearing in the linearized col- 
lisional term are, respectively, the absorption and scattering cross sections in D 
spatial dimensions, N is the number of scatters and we have considered coherent 
scattering only. We do not evaluate them here, since for our purposes this is not 
relevant. The phase function p(n,n') is given by relation 



da s (n.n) = p(n, n')a s dflu . 



(74) 



If we consider only Thomson scattering, (a s = cr Th ), we have in D-dimensions 
(see Appendix B) 

1 



p(n, n') 



S r 



1 + , ° ,J n"n v + —br)M u ri v + —hau) 
(£) — l) 2 D D 



(75) 



which reduces to the usual Thomson formula for D = 3 (see |p4,25|). 

In order to compute the integral in (73), we consider the total distribution func- 
tion as given by (1), (2), and (5), namely 



F = F (0) + A + BX + C^inV + ^W v ) . 
Inserting (75) and (76) into (73), and evaluating the integral, we get 



(76) 



k^F = -u {Na a (u)(l - e~ uj/T )A + N [a a (u)(l - e^ /T ) + 



+N 



, n/, u/t-, D(D 2 -3) 



(77) 
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Comparing the above expression with equation (8) we obtain the transport co- 
efficients 

K (uj) = Na a (uj)(l-e-" /T ) , (78) 
K 1 {u) = N[a a {u){l-e-"l T ) + a Th ] , (79) 

, 2 ( w ) = iVKH(i- e -/V^^sJ • (so) 

Note that when the Thomson scattering dominates, kq(u) is negligible, thereby 
implying that there is no radiative bulk viscosity in this case (see (25) and (41)). 
Under such a condition an homogeneous and isotropic mixture of matter and radi- 
ation could expand in equilibrium. However, for an inhomogeneous and anisotropic 
medium, heat conduction and shear viscosity are always present . The Rosseland 
means for K\ and ki reads 

k = »k • (81) 

1 1 (D- l) 2 (D + 2) 



82) 



k 2 N D(D-3)a Th 

It thus follows that the transport coefficients (41), (42), and (43) are given by: 

D + 1 a„ mn 
Xd = „ AT ° T D , (83) 



D No. 



Th 



_ (D + l)(D-l?a D T^ 

Id D 2 (D 2 - 3) Na Th ' 1 1 



C=0 • (85) 

In particular, for D = 3, the coefficients first computed by Straumann || are recov- 
ered. 
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VI. CONCLUSION 



As we have seen, kinetic theory is indeed a very powerful approach to describe 
radiative plasmas irrespective of the number of spatial dimensions. It provides the 
transport coefficients and the standard entropy production rate, which is obtained 
from non-equilibrium thermodynamics. In addition, through the same algorithm 
used to define the Eckart temperature, the main macroscopic quantities are also 
consistently settled to have the same form of equilibrium, as required by the first 
order nonequilibrium thermodynamics. 

The Eckart approach considered here is plagued with serious undesirable features, 
like instabilities of the equilibrium states and superluminal velocities of the thermal 
and viscous signals [p6fl - [pSfl . Such shortcomings can be circumvented by the so- 
called second order thermodynamic theories as developed by Muller [E(J and Israel 



27]. In this way, it is interesting to reconsider the kinetic treatment developed 



here in the framework of the causal or transient relativistic kinetic theory . The 
corrections of the extended theories, however, will be important only in regimes 
where the mean free paths are comparable to the macroscopic length scales |29| . 



Finally, we remark that is not difficult to widen the scope of the results de- 
rived here to include general relativistic effects. As usual, one must assume that 
the mean free time of each process is much smaller than any typical scale of the 
gravitational field, and at the level of the computations, to implement the minimal 
coupling, viz.; to replace the ordinary derivatives by covariant derivatives and the 
(D+l)-Minkowski matrix 7]^ by g^ v ||. Naturally, D-gravity does not change the 
form of the D-dimensional transport coefficients, which should be considered in the 
analysis of the entropy production in high dimensional cosmologies, as well as in 
the role played by dissipative processes on the problem of dimensional reduction. In 
particular, as shown in section 3, the radiative bulk viscosity does not play any role 
in a D-dimensional radiation-dominated FRW universe. This generalizes the well 
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known results obtained by Weinberg in the three dimensional case 0. Qualitatively, 
significant bulk viscosity may take place only if the plasma does a transition, for in- 
stance, from radiation to a matter-dominated phase. However, if only the Thomson 
scaterring is considered, the bulk viscosity also vanishes, regardless of the number 
of spatial dimensions. Further considerations in these topics will be presented in a 
forthcoming communication. 
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APPENDIX A: IRREDUCIBLE POLYNOMIALS NORMALIZATION 

For completeness, we now outline the proofs of formulas (11)-(13), which are the 
key ingredients to perform integrals defining the usual thermodynamic quantities. 
These formulae on the normalization of the irreducible polynomials are frequently 
cited in the case D = 3 (see, for instance, ref(s) [4], [5] and [7]). Nevertheless, as 
far as we know, there is no simple proof of them in the general case. As mentioned 
earlier (see section 2) we make all calculations in the comoving frame. 

The hypersphere parametrization is given by 

k 2 + k\ H h k 2 D = uj 2 

where 
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uj cos 6 



i 



uj sin 61 cos 6 2 



uj sin 61 sin 6> 2 cos #3 



k D -i 



uj sin ^! sin 6 2 . . . sin 6 



D-2 



cos 6 



'D-1 



k D 



uj sin 61 sin 9 2 ■ ■ ■ sin 6 



D-2 



sin 6 



D-1 



with < 9j < it for j = 1, 2, . . . , D — 2 and < 6d-i < 27r. The element of volume 
on S D1 is given by 

dVtu = sin D - 2 9 1 sin D ~ 3 6 2 . . . sin 1 6 D _ 2 sin B D - X dQ x dQ2 ■ ■ ■ d& D - X , 
and the volume of the (D — l)-sphere is given by equation (10). 

a) Proof of equation (11) 

In the comoving frame uj = k° and n° = 0, thereby implying that (11) is valid 
for ji — and v = 0, 1,2,3 (both sides are identically zero). Now consider % ^ j. 
Since in this case h l i = 0, we need only to prove that f s ^ n % n?d£lu = 0, or the 
equivalent formula J s ^ k' l k j dilu = (note that k % = uon 1 ). Assuming, without lost 
of generality, that i < j, and using the above parametrization of the (D — l)-sphere, 
it follows that 



where A is a constant which appears from a multiple integral involving all angular 
variables but Q; t . 

For i = j we have from spherical symmetry 




(Al) 




(A2) 
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and since h n — — 1, % — 1, 2, . . . , D we can write 

/ n *dSl u = -±h«S D _ 1 , (A3) 
which means that (11) is valid for /i, v — 0, 1, 2, 3. 

Proof of equation (12) 

By using equation (11) we get from the left hand side of (12) 
— *— / (n p n u + l/i"") {rfn p + ±-h ap )dn v = 

■nf i n v n a n"dn u - -^-h pv h ap . (A4) 



^ S D-1 D 2 

Let us now consider the integral 

— !— f n p n u n CT n p dtt u . (A5) 
j d _ 1 Js D1 

To evaluate it we observe that h pu is the only second rank D-dimensional tensor at 
our disposal which does not depends on the vector n p , and that all permutations of 
(fj,uap) leave us with the same expression. The only 4-rank tensor matching these 
requirements is given by h pu h ap + h pa h vp + h pp h ua . The integral must be proportional 
to this tensor, and the proportionality constant should depend only on the dimension 
D. Therefore we can write 

- — f n lt n v n <T n p dn u = f(D)(h^h ap + hTh vp + h pp h ua ) , (A6) 

^D-l ^ S D-1 

where f(D) is the constant to be determined. Since we have tensorial equation, we 
can evaluate the scalar f(D) in the comoving frame. By specializing the integral for 
n p = n u = rf = n p = n 1 = cos#i, we obtain 

f(D) = — L_ / nSdttu . (A7) 

OD D _ 1 JS D _ 1 

A straightforward (although a bit lengthy) calculation leads to 

!{D) -D(h¥ } • (A8) 
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and using (A. 8), (A. 6) and (A. 4) one immediately obtains (12). 



Proof of equation (13) 

In the comoving frame (n° = 0), we need to proof (13) only for the spatial 
components. Suppose that we have a odd number, say p, of components of the 
vector n. Since rt = k l /k° we get 

f n l n j . . . n l n r dQu = f k l k j . . . k l k r dQ v . (A9) 

The vectors k = (k°, k 1 , k 2 , . . . , k D ) and k = (k°, —k 1 , —k 2 , . . . , —k D ) cover the same 
D-dimensional sphere S D _ 1 . It thus follows that 

f k*k j ...k r dQu = {-l) p f k { k? ...k r dQ v . (AlO) 

Since p is odd we get J s ^ k' l k j . . . k l k r d£lu = and so the proof is completed. 

APPENDIX B: D-DIMENSIONAL THOMSON CROSS SECTION 

In order to evaluate the phase function p(n, n') one must impose some hypothesis 
reflecting the symmetry of the specific scattering problem. For Thomson scattering, 
the results of references [5] and [19] are easily generalized as follows: 

HI p(n,n') is written in terms of irreducible polynomials. 

H2 p(n,n') — p(n' ,n) 

These two hypothesis imply that p(n, n') is constructed out from products of an 
even number of irreducible polynomials. Thus, we can write 

pin, n') = A + Bn a n' a + C{n»n v + ^hT) {n'^n' u + -^) , (Bl) 

where A, B, and C are constants to be determined. 
H3 p{n,n') is normalized such that: 

/ p[n,n')dSl v = \ . (B2) 

^D-l 
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H4 The mean value of n M on the sphere S D _ 1 is zero, that is: 

f p(n,n)n'^d% = . (B3) 

Substituting (B.l) into (B.2) we get A = — ^ — , and from (B.3) we get B = 0. Hence, 
p(n, n') reduces to 

p(n, »') = ^ + + (« + ^ V) • (B4) 

To compute the coefficient C, we consider the comoving frame where n = (0, n) = 
(0, rii, 1^2, ■■■ , Tin), n' = (0, n') = (0, n[, n' 2 , . . . , n' D ), and define the angle 9 by 

n^n'^ = -n t v! i = cos 9 . (B5) 

Inserting (B.5) into (B.4) we get 

p(n,n') = -^—[l + a(cos 2 9-h} , (B6) 

£>— 1 U 

where a = CS D _ 1 . The differential scattering cross section is defined by (see Fig 1) 

da = sin 2 (ftdflu = p(n, n')a s dVLjj , (B7) 

while the total scattering cross section a s reads 

r D - 1 

<J S = I da = —fiS^ . (B8) 
J s n ., u 



J D-\ 



It thus follows from (B.6), (B.7), and (B.8) that 

sin 2 0= ^-i[l + a(cos 2 ^--^)] . (B9) 

From Fig. 1 we get the relation between the angles ( see, for example, reference [20] 
for a 3- dimensional version of our argument) 

cos <fi = sin 9 cos ip . (B10) 

Now, consider the sphere S D _ 2 C S D _ 1 , such that n _L S D _ 2 , that is, n is orthogonal 
to all vectors describing S D 2 . By denoting < A > as the mean on S D 2 of an 
arbitrary quantity A, namely 
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<A>=- 5 —l AdVL D _ 2 , (Bll) 



Q Q 

w D-2 °D-2 

we obtain (note that only the angle ip relies on S D 2 ) 

1 



< cos 2 i) >= j^—^ , (B12) 
and from (B.10) and (B.12) 

< sin 2 >= 1 - - sin 2 9 . (B13) 
Now, taking the average on S D _ 2 of equation (B.9) and using (B.13) we get 

c =irrw^sz • (B14) 

and inserting (B.14) into (B.4) we find the D-dimensional Thomson cross section 

P(n, n>) = + D _ (nV + ^r)(n;< + . (B15) 

As expected, for D = 3 this expression reproduces the usual Thomson formula 



PTh (n,n , ) = -^(l + co S 2 9) . (B16) 

107T 
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FIG. 1. Thomson Scattering in D dimensions. The geometry of the scattering is 
the same as in 3 dimensions but the mean of the angle tp has been done on the 
(D — 2)-dimensional sphere S D 2 . 
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